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I point out how coherence vortices, i.e., topological defects in a correlation function, could help 
explore new physics if they are created in matter waves. Vortex dynamics could be studied in up 
to six dimensions, and spin topological defects unseen in lower dimensions could be created. A 
rudimentary proof-of-principle experiment is sketched and simulated, in which three Bose-Einstein 
condensates are used to create and detect coherence vortices. 



Coherence vortices is a concept first coined in optics 
research [l], [1| . These are phase singularities that exist 
not directly in a physical field, but in a spatial correlation 
function of a field. Typically, one may have a partially 
coherent field, such as a laser beam with limited coher- 
ence time, where the observed expectation value does not 
exhibit a phase singularity, but a double-slit interference 
experiment may reveal a phase winding around a sin- 
gularity in the field-field correlation function. As such, 
coherence vortices are defined in a coordinate space with 
twice the dimensionality of the physical space. Simula 
and Paganin pointed out that this means that coherence 
vortices are possible in a one-dimensional system [3] ; fur- 
thermore, they alluded to the possibility of observing co- 
herence vortices in matter as well as in optics. 

A vortex in a two-point correlation function suggests 
the intriguing prospect of studying vortex dynamics in up 
to six dimensions. However, for the study of dynamics, 
optical systems are impractical and one needs to turn to 
coherent matter waves - partially or fully Bose-Einstein 
condensed gases of ultracold atoms. Matter waves have 
several advantages as well as disadvantages with respect 
to optics when it comes to studying coherence vortices. 
Detection is a very different affair in matter waves, but it 
can be done, as will be discussed in detail below. Among 
the advantages are not only the possibility to study dy- 
namics, but also the prospect of using multi-component 
condensates for creating higher-dimensional topological 
objects. A four- or six-dimensional space will offer new 
possibilities for these. 

The purpose of this paper is to point out the new 
prospects for research that are made possible by co- 
herence vortices in matter systems; notably, high- 
dimensional vortex dynamics and high-dimensional topo- 
logical defects. In addition, I sketch and simulate a 
clumsy but doable scheme for the realization and detec- 
tion of coherence vortices in a system of trapped cold 
bosons. I see this as a proof of principle, which helps 
open up the way towards six-dimensional vortex dynam- 
ics. 

Let us first note that there exist already theoretical 
proposals that would in effect contain coherence vortices; 
namely the quantum fluctuating vortices studied by Mar- 



tikainen and Stoof [4] , and the slowly-rotating small Bose 
and Fermi systems studied by Reimann, Kavoulakis et al. 
{H EJ . What these systems have in common is that the 
exact ground state would contain a vortex whose posi- 
tion is subject to quantum uncertainty. The symmetry 
will be broken and a definite vortex position is chosen by 
any finite symmetry-breaking potential; this destroys the 
superposition and thus the coherence vortex is reduced 
to an ordinary vortex. The scheme I describe below is 
more robust, though perhaps less elegant. 

System. - Cold bosonic atoms are described by com- 
muting field operators ^a( r ) 7 V4( r )> where a designates 
species and/or spin state (F, mp), and r is the spatial co- 
ordinate vector. If the atoms form a Bose-Einstein con- 
densate (BEC), the expectation of a single field is non- 
vanishing, ^ a ( r ) = (^o(r)), and its expectation value is 
called the condensate wave function. The vortices sup- 
ported by such wavefunctions are, of course, well studied 
[7j. Instead, we are here interested in phase singularities 
in the field-field correlation function, 

5oo (r,r') = (^(r)^ (r')>, (1) 

which is a function of 2 x D coordinates, where D is the 
number of dimensions in which the atoms can move. Such 
phase singularities may exist also when there are none in 
\I/ a (r). Henceforth, we will be dropping the subscript a 
when we are dealing with only one component. 

Topological defects. - In a single-component fluid in D 
dimensions, the only type of topological defect that can 
exist is a vortex, characterized by an integer phase wind- 
ing number around a closed loop. Such a loop surrounds 
a phase singularity of dimension D — 2, i.e., a point vortex 
in two dimensions or a vortex line in three dimensions. 
In D — 4 the singularity lies on a two-dimensional (2D) 
surface, and it is helpful to imagine the full 4D space 
as a succession of 3D spaces, each (except in degenerate 
cases) containing a vortex line, whose trajectory varies 
continuously as we move between adjacent 3D spaces. A 
closed loop is drawn around the vortex line in one of the 
3D spaces. If the loop is displaced by a sufficiently small 
amount in any direction, including the fourth, it will still 
encircle the singularity. Thus the 2D sheet retains the 
topological stability we associate with a vortex. 
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In a multi-component quantum fluid, higher- 
dimensional topological defects are possible, as is 
well known in the 3D case. The allowed topological 
defects in a given quantum fluid are classified using 
homotopy theory (j|-[lQ|: Let G be the symmetry group 
of the order parameter, and let H be the group of 
operations that leave a specific ground state invariant. 
Then the order parameter space M = G/H is is one-to- 
one correspondence with the set of distinct degenerate 
ground states of the system. Now the first homotopy 
group of M, 7Ti(M), gives the allowed topological defects 
that may be encircled by a ring; these are (singular) 
point defects in 2D, line defects in 3D, and planes in 4D. 
Likewise, the second homotopy group 7r2(M) enumerates 
those topological defects that can be encircled by a 
spherical shell, and in 3D those are point defects, i.e., 
monopoles. For instance, in a single-component BEC, 
M = S 1 , the circle; and we have tti(S 1 ) = Z, where Z 
is the set of integers. These integers are the possible 
winding numbers for a vortex. It is now plain to see that 
in higher dimensions, there may be topological defects 
classified by higher-order homotopy groups, which do not 
make sense in lower dimensions - e.g., its characterizes 
singular point defects in 4D or sheet defects in 6D. In 
a single-component BEC, again, all homotopy groups 
higher than 1 are zero; ^(S 1 ) = tts(S 1 ) = etc; which 
implies that the only topologically stable defect for such 
a fluid is the vortex. For the ferromagnetic state of a 
spin-1 condensate, M=SO(3), and its homotopy groups 
are tt 3 (M) Z and tt 4 (M) tt 5 (M) Z 2 , where Z 2 
is the set of integers modulo 2. Further investigation 
of the homotopy groups associated with the correlation 
function matrices g a b of the ground states in spinor 
condensates must be left to future studies. 

Creation. - We now turn to the experimental creation 
and detection of coherence vortices. A number of ingre- 
dients are necessary for the creation, and a few more for 
the detection, of coherence vortices. First, one needs a 
non-trivial correlation function. If the whole system is 
fully coherent - i.e., one pure BEC - vortices exist in 
the correlation function if and only if they are there in 
the condensate wavefunction itself, and the study of the 
associated coherence vortices would be pointless. A par- 
tially coherent cloud, such as bosons in a one-dimensional 
or disordered potential, will also not present coherence 
vortices in most cases, but only a decaying correlation 
function. Instead, I sketch an experiment that uses the 
random relative phase of two independent BECs. 

The experiment starts out with a pair of independent 
BECs, for simplicity in two dimensions. Pairs of indepen- 
dent BECs are created by holding atoms in two separate 
traps while cooling them down below the critical temper- 
ature [11]. Such pairs of BECs will produce interference 
patterns when brought together, but the two are statis- 
tically independent in the sense that their global phase 
difference - and thus the phase of the interference pat- 



tern - varies randomly from shot to shot. Interference 
between the two independent BECs is, thus, washed out 
only in the average over many experimental runs. In 
an optical system, a beam has a finite coherence time 
and thus it self- averages. It is possible that one could 
construct something similar by making clever use of the 
finite coherence length of an elongated or disordered sys- 
tem, but I will not follow that line of thought here. I 
thus opt for averaging over repeated single shots as the 
straightforward (but tedious) way to construct an en- 
semble average - by hand, so to speak - in this proposed 
experiment. 

The BECs are put in different spin states and released 
towards each other head on. The two spin states a will be 
labeled 1 and 2, respectively, and the trapping potential 
for each state is written as 



V a (r) = ^mu 2 (r - r 0o ), 



(2) 



where the trap center i*o a is taken to be different for the 
two BECs. They are first Bose-Einstein condensed at 
a distance large enough to avoid contact, and then adia- 
batically brought to their starting positions, ro a - At time 
t = the traps are suddenly moved to target positions 
r' Qa . The BECs are accelerated towards these new poten- 
tial minima, and the resulting collision excites vorticity. 
At a final time t, the atoms of species 1 are transferred 
into species 2, and the resulting state is recorded. Di- 
rect observation of the density profile will reveal a set of 
vortices, whose positions depend on the initial, random 
phase difference between the two BECs. We denote this 
phase difference 0. It is in the transfer at time t that 
this phase difference will be decisive. 

The process is simulated for a binary mixture of 4.4 x 
10 4 atoms of 87 Rb equally distributed among the states 
\F, Trip) = |1, 1) and \F, mp) — |1, —1), each in a trap of 
frequency uo = 2tt • 20 s _1 in the plane, and in addition 
a tight trap of frequency 2tt • 400 s _1 is assumed along 
the z direction; this direction is integrated out so that 
the physics in the 2D plane is simulated. The starting 
positions are at roi = (12 /im, 0), and 1*02 = —1*01, and 
the target positions are both at the origin, = r f 02 = 0. 
This system is described by the two-component Gross- 
Pitaevskii equation, 
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with U a b = 47r?i a a &/m, where the s- wave scattering 
lengths are an = a22 = 100. 4ae and a\i = 101. 3aB, 
where an is the Bohr radius. A definite value of 6 is 
used in each simulation run. Two examples, for different 
O, are shown in Fig. [TJ The physics of the violent colli- 
sion is seen to be mainly governed by pass-through [12] 
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FIG. 1. [Color online] Collision of two BECs for two cases with 
different and randomly chosen initial phase difference. Simu- 
lation parameters as in text. Density and phase are recorded 
at t = 100 ms. (a), (c): Density; (b), (d): Phase, (a), (b): 
6 = 1.98296 radians; (c), (d): G = 4.19545 radians. Units 
for the density is such that its maximum value is 1. Posi- 
tions of vortices are indicated with crosses in the figures, and 
ant ivort ices with rings. 

and quantum interpenetration [l3[. Subsequent trans- 
fer of the population to component 1 results in a single 
BEC containing a set of vortex- ant ivortex pairs whose 
positions - in the Gross-Pitaevskii picture - depend on 
6. Averaging over the phase difference will wash out 
the vortices. However, the vorticity in the correlation 
function will not. The exact correlation function after 
transfer is simply 

g(r,r>) = tf? (iWr') + ^(r)* 2 (r'), (5) 

since the cross-term vanishes for independent conden- 
sates. The vortices in g(x, y, x' ', y') corresponding to Fig. 
[T]are imaged in Fig. [2] as a sequence of 3D frames in the 
space (x, y, x'), with a fixed y' for each frame. There are 
vertical vortex lines visible in all four isosur faces; these 
are inherited from the vorticity created in the collision 
(cf. Q). In addition, a tangle is created of vortices not 
perpendicular to the xy plane; it is tempting to think of 
these as more genuine coherence vortices. 

Note that before the transfer, the vorticity in the two- 
component system takes the form of skyrmions or filled- 
core vortices. However, it would be incorrect to call them 
coherence vortices, since their locations are independent 
of the phase difference 6 and thus they are manifest di- 
rectly in the density itself, even after averaging. 

Detection. - In order to observe the coherence vor- 
tices experimentally, one must reconstruct the field-field 
correlation function g(r, r'). For cold atoms, there exists 
no practical method for measuring it at arbitrary coor- 
dinates. For optical systems, one may let the light pass 
through two pinholes and observe the phase of the re- 
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FIG. 2. Vortices (isosurfaces of the phase winding) of the 
correlation function g(x,y 1 x / ,y'), at three-dimensional slices 
of the four-dimensional space, fixed at (a) y — 0, (b) y — 
3.42/xm, (c) y' — 6.83/xm, and (d) y' — 10.25, respectively. 
The parameters are as in Fig. [1] 

suiting interference pattern. Ref. [3] mentions this type 
of method as an alternative also for matter systems, but 
does not discuss how. It is possible that a scheme could 
be devised where atoms in small neighborhoods of two 
points are excited and then brought together. Instead, 
one may obtain an approximate correlation function by 
experimentally imitating the type of averaging that is 
done in c-field calculations [lij . 

In a single realization, the spatial profile of the con- 
densate's phase can be measured by letting the BEC in- 
terfere with another one initially held at a distance. (In 
the present scheme, this means managing three BECs, 
whereof two should be independent.) One can now record 
the set of phase differences between all points in the sin- 
gle run, and then average this function over repeated runs 
to obtain the approximate correlation function, 

g^ p (r,r') = (e^-^), (6) 

where 0(r) is the phase at point r inferred from the fringe 
pattern in a single experimental run. The approximation 
built into this procedure is the neglect of spatial density 
variations, as seen in Eq. (j6]). It is often a good one as far 
as the locations of the centermost vortices are concerned. 
We simulate an actual experiment by making 30 individ- 
ual runs with random initial values of and random 
noise on the wavefunctions \£j on a 5% level. In Fig. [3j 
the simulated correlation function g exp is compared with 
the exact one that includes density information. The re- 
sult may be improved by reading off a coarse-grained den- 
sity profile from the interference fringes, if one so wishes. 

The time development of the coherence vortices may 
now be studied by doing repeated experiments; since we 
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FIG. 3. [Color online] Field-field correlation function g(r, 1*0), 
with 1*0 fixed. The position of 1*0 is indicated with a star, (a), 
(b): Exact correlation function, (c), (d): Approximate cor- 
relation function of Eq. (0, which simulates the result of the 
proposed experiment, (b), (d): Phase of g; (b), (d): Absolute 
value \g\. Normalization is such that the maximum value of 
\g\ is 1. The parameters are as in Fig. [U 



are studying average quantities, the destructive imaging 
technique is not a problem. 



FIG. 4. [Color online] Time development of the field-field 
correlation function g(r, r ), with 1*0 fixed. The modulus 
|g(r, r )| is displayed, (a): Time t = 25ms, (b): t — 50ms, 
(c): t — 75ms, (d): t = 100ms. Positions of coherence vor- 
tices are indicated with crosses in the figures, and ant ivort ices 
with rings. The position of 1*0 is indicated with a star. The 
parameters are as in Fig. [T] 

tion function mirrors the physical process: two conden- 
sates approaching, colliding and interpenetrating, creat- 
ing vortex- antivortex pairs in the process. This picture is 
made especially clear by choosing the 2D plane of view to 
be parallel to the physical space; as we have seen, vortex 
tangles move in all four dimensions. 

Concluding remarks. - Summing up, this Letter wishes 
to point out two (related) reasons for studying coherence 
vortices in cold matter systems: It opens up for study- 
ing dynamics of vortices in high dimensions, and new 
types of high-dimensional topological defect, unseen in 
3D, can be created here. As a specific example, a system 
of two colliding BECs was simulated, in which a tangle 
of two-dimensional coherence vortex sheets was created. 
A recipe for a proof-of-principle experiment was laid out 
where coherence vortices are created and detected using 
current experimental technology. 

The author is indebted to Harri Makela for help with 
homotopy groups and to Jani Martikainen for feedback 
on the manuscript. 



Correlation functions for four time instances are shown 
in Fig. [4j It can be seen how the behavior of the correla- 
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